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Results are presented from a time— dependent, numerical investigation of super- 
Eddington spherical accretion onto black holes with different initial conditions. 
We have studied the stability of stationary solutions, the non-linear evolution of 
shocked models and the time-dependent accretion from an expanding medium. 



1 Introduction 

Stationary spherical accretion onto black holes is certainly a well-known and 
extensiveli^ studied topic (for a comprehensive review see Nobili, TuroUa and 
Zampieri B, Zampieri, Miller and TuroUa □). Up to now all time-dependent 
investigations of spherical accretion, have been concerned either with_the anal- 
ysis of the stability of isothermal B, isentropic □ and optically thiclO' □ flows 
or with the definitioiLol the parameter space within which high temperatuxe 
solutions might exist Q' El. Very recently, Colpi, Shapiro and WassermanEj 
have studied time-dependent hypercritical accretion of an initially expanding 
medium to describe the late fall-back of material onto the compact remnant 
after supernova explosion. 

Despite the fact that spherical accretion has been extensively investigated, 
mainly for shedding light on the efficiency of the radiation generation, sev- 
eral aspects require further consideration such as, for instance, investigating 
the stability properties of high temperature solutions and searching for the 
existence and non-linear evolution of possible heated or shocked models. The 
results of this calculation will be presented in Sec. 2. This study will serve 
also as a starting point for the investigation of very interesting astrophysical 
problems, such as the late fall-back of material onto a supernova remnant. A 
preliminary analysis of this problem will be presented in Sec. 3. 



2 Stability of stationary solutions and evolution of shocked models 

We have studied the evolution of a self-gravitating, perfect hydrogen gas in 
the gravitational field of a non-magnetized, non-rotating black hole within the 
framework of time-dependent radiation hydrodynamics. The present investi- 
gation has been carried out by means of a general relativistic treatment of the 
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Figure 1: The gas temperature T and the radiative luminosity I (in Eddington units) are 
plotted versus log(r/rg) (rg is the gravitational radius) at different times for a HL model 
with initial m = 28.3 (in Eddington units). 



radiative transfer equation, which exploits the expansion of the specific inten- 
sity of the radiation field into moments til. A detailed review of the derivatiori 
of the fundamental equations is presented in Zampieri, Miller and TuroUaa. 
The complete system of time-dependent, radiation hydrodynamics equations 
plus the Einstein Field Equations has been solved using_an original numerical 
method, based on a Lagrangian finite difference scheme 

Depending on the physical conditions in the accretion flow, stationary so- 
lutions of hvpercritical spherical accretion onto black holes show very different 
behaviours □. In the LL (Low Luminosity) models, the inner region of the 
accretion flow is in LTE with radiation at a low temperature (T ~ f ) 
and luminosity and efficiency turn out to be very small (L < 4 x 10'^^ erg s~^, 
e ^ fO"''). The HL (High Luminosity) models are characterized by very high 
inner temperatures [T ^ 10^-fO^° K). In the intermediate region (fO^-fO^rg, 
where rg is the gravitational radius), Compton cooling and free-free emission 
nearly balance Compton heating, producing a total luminosity 2-3 orders of 
magnitude larger than that emitted by LL models with the same density. 

The time-dependent analysis of LL models shows that they are stable to 
both thermal and radiative perturbations □. Much more spectacular is the 
phenomenology shown by HL models □. At high accretion rates (m ^ fO, m 
in Eddington units), a thermal instability appears around 10'^ after about 
2-3 s, as can be seen from Fig. f, and the temperature increases by almost 
an order of magnitude. This instability is due to the fact that the radiative 
heating rate is greater than the radiative cooling rate and, at the same time, the 
heating time is shorter than the dynamical time. A few seconds after this, the 
velocity profile starts to deviate significantly from free-fall owing to the large 
drag exerted by the internal pressure gradients. A compression wave develops, 
which becomes progressively steeper as it propagates outward and, after 8-10 s. 
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Figure 2: The gas velocity v (in units of c^) and the gas density p (in units of po) are 
plotted versus logi? (R = r/r^, rJJ is the initial accretion radius) at different times. Initial 
data refer to SN 1987A; po ~ IQ-^ g cm'^, c° = 350 Km s-\ = 1.6 X lO" cm. 



a hydrodynamic shock forms. The luminosity profile has the typical behaviour 
shown in Fig. f with the long-dashed line, showing a significant initial transient 
increase lasting 8 s. The shock front moves outward at an approximate speed 
of 10* cm/s ~ lO^^c. Hence this solution is definitely non-stationary. 

HL models with low accretion rates {ifi < 10) do not show any thermal 
instability on comparable evolutionary timescales. Having lower densities, the 
radiative heating and cooling are comparatively less efficient than compres- 
sional heating and the gas is essentially adiabatic up to large radii. 

3 Spherical accretion from a uniformly expanding medium 

After (Type II) supernova explosion of a high mass star, a reverse shock may 
arise from the collision of the outflowing gas with a sufficiently massive outer 
stellar envelope, leading to the formation of a dense, expanding flow sorround- 
ing the compact remnant. Models for SN 1987A show that the core density 
is roughly uniform {p ^ 10~^ g cm~^) and that the initial velocity v oc r 
with a maximum value ~ 2000 Km s'^. After the expansion phase, gravita- 
tionally bound shells can accrete onto-jthe compact remnant at hypercritical 
rates. Colpi, Shapiro and Wassermanll3 found that the total accreted mass is 
^ O.IM0 and could suffice to drive a remnant borne as a neutron star to further 
collapse to a black hole for a relatively soft equation of state. In this respect 
we note that, for the large accretion rates occurring during the early phase, 
accretion onto a neutron star might not be drastically different from accretion 
onto a black hole because of the trapping of radiation and the jeffectiveness of 
neutrino losses close to the neutron star surface (see Chevalierli3) . |_ 
Starting from the initial conditions and parameters expected for SN 1987AI1: 
we have repeated the calculation of Colpi, Shapiro and Wasserman for a spheri- 
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cally accreting polytropic gas with T = 4/3, including self-gravity. The numer- 
ical results are shown in Fig. 2. As can be seen looking at the velocity profile, 
at the beginning of the evolution (dash-dotted line) three distinct regimes 
can be recognized: free-fall in the region below the accretion radius r^, near 
hydrostatic equilibrium between ra and the marginally bound radius rmb (be- 
low which mass elements are gravitationally bound) and uniform expansion 
above rmb- At t ~ 4700 s, the accretion rate reaches its maximum value, 
m ~ 2 X 10^". The late fall back can be described in terms C|£_a dust infall 
in which only those shells within r^b are eventually accreted lid. The total 
accreted mass is Alacc — O.GSA/q, in agreement with the value estimated by 
Colpi, Shapiro and Wasserman. Then, although the initial mass contained in 
the expanding shells is ~ 5Mq, only a small fraction is accreted. Since for the 
Birkhoff theorem the mass outside Tq has no effect on the motion of the inner 
gas and since Mace is very small in comparison with the mass of the central 
black hole, self-gravity does not play any important role. 

The present simple hydrodynamic model of spherical accretion from a uni- 
formly expanding polytropic gas represents the starting point of this investi- 
gation and does not include the transfer of radiation, the role of which might 
probably become important. In fact, although matter and radiation are pre- 
sumably in LTE and the efficiency of accretion onto black holes is very low, the 
actual value of the accretion rate is so high that near Eddington luminosities 
might be produced, strongly influencing the dynamics of the accretion flow and 
hence the total accreted mass. This study along with a companion analysis of 
spherical hypercritical accretion onto neutron stars are presently under way. 
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